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1. Introduction 

A new generation of experiments devoted to hadron spectroscopy, including GlueX at Jef- 
ferson Lab, PANDA at GSI/FAIR and BES III intend to make measurements with unprecedented 
precision and in previously unexplored mass ranges and quantum numbers. Both in meson and 
baryon spectroscopy, there are many experimentally observed excited states whose physical prop- 
erties are poorly understood and could use theoretical input from lattice QCD to solidify their 
identification. Aside from masses, other excited-state quantities that could be computed on the 
lattice, such as form factors and coupling constants, would be useful to groups such as the Excited 
Baryon Analysis Center (EBAC) at Jefferson Lab, where dynamical reaction models have been de- 
veloped to interpret experimentally observed properties of excited nucleons in terms of QCD [jl|]. In 
certain cases, input from the lattice may be helpful in determining the composition of controversial 
states, which may be interpreted as ordinary hadrons, tetra- or pentaquarks, hadronic molecules or 
unbound resonances. 

Among the excited nucleon states, the nature of the Roper resonance, Af(1440) P\\, has been 
the subject of interest since its discovery in the 1960s. It is quite surprising that the rest energy of the 
first excited state of the nucleon is less than the ground-state energy of the nucleon's negative-parity 
partner, the Af(1535) S\\ [0], a phenomenon never observed in meson systems. There are several 
interpretations of the Roper state, for example, as the hybrid state that couples predominantly to 
QCD currents with some gluonic contribution or as a five-quark (meson-baryon) state [J3j] 

Due to the greater impact of systematic errors, such as finite-volume effects and discretization 
errors, on excited states, there is an essential need to examine calculations of excited masses more 
carefully during analysis. It is common in such calculations to use the ground-state masses, such as 
nucleon mass, as a starting point for analysis but without further checks of the consistency of the ap- 
proach. The nucleon mass has been demonstrated with good consistency (a great demonstration of 
the universality of lattice QCD using different fermion actions by various groups with independent 
analyses) in both quenched and dynamical Nf ■ = 2 + 1 cases. For dynamical ensembles, there are 
more variables in terms of algorithm, scale setting, etc., but an approximate universality is achieved 
among different groups (see middle panel of Fig. [[]). However, beyond the ground state, there ex- 
ists a diverse distribution of excited-state masses as functions of m n . There are big discrepancies 
in the calculated nucleon first-excited mass, creating an apparently chaotic atmosphere, as shown 
in the left panel of Fig. [j]. Note that the errorbars here are just statistical, none of the systematic 
errors (such as quenching, finite-volume effect, etc.) are estimated. We re-address the same issue 
by modifying the axes to be in terms of the dimensionless quantity ML (as shown in right panel 
of Fig. |]); we find that the Roper masses are roughly inversely proportional to lattice size. Now 
we see a better agreement (or universality) among the lattice QCD Roper-mass calculations; most 
of the Roper masses agree within 2 standard deviations of the numbers in Ref. [ft]. This suggests 
that finite-volume effects can be more severe for excited states than ground states and that careful 
examination of such systematic errors is crucial. Similarly, we cannot ignore other systematic er- 
rors that may arise even if we did not observe the effect in the ground state. Excited-state analyses 
should proceed with greater caution. 

To tackle the challenge of extracting reliable excited-state energies, the Hadron Spectrum Col- 
laboration (HSC) has been devoting resources and effort into resolving the mystery from the funda- 



2 



Excited-Nucleon Spectroscopy with 2+1 Fermion Flavors 



Huey-Wen Lin 



* 1.5 



*{>%*> 

tf* t * I 



; ! ' 



It* 

* Mathur (3.2 fm) 
+ Sasaki (3.0 fin) 
o Lasscock (2.6 fm) 
a Burch (2.4 fm) 
□ Brommel (2.4 fm) 
Basak (2.3 fm) 



A# (GeV- 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 



■ HSC 

A RBC-UKQCD 

: LHP 

o PACS-CS 
O MILC 




+ Mathur (3.2 fm) 

+ Sasaki (3.0 fm) 

o Lasscock (2.6 fm) 

- Burch (2.4 fm) 

□ Brommel (2.4 fm) 

< Basak (2.3 fm) 



Figure 1: Summary of published quenched (Nf = 0) lattice QCD calculations of the nucleon and Roper 
masses in GeV (left) and in terms of the dimensionless product of the Roper mass and lattice size L (right). 
The middle panel shows the dimensionless products of lattice nucleon mass and rrj as a function of m^ro 
from Nf = 2+1 ensemble. 



mental point of view. It is difficult to extract the excited state reliably without sufficient information 
carried in the nucleon correlators. We need correlators which clearly contain the masses of specific 
quantum numbers and allow us access to data points that reach higher excited states before the 
signal decays exponentially away. To successfully and reliably extract these excited-state energies, 
we need better resolution in temporal direction (large superfine isotropic lattices or anisotropic lat- 
tices), separation of the signals for individual states (variational method), and most importantly, 
operators that have good overlap with various desired quantum numbers (cubic-group irreducible 
representations and some way to put in many operators that contribute unique signals). 

2. Nf = 2 + 1 Clover Anisotropic Lattices and Updates 

To improve our ability to extract higher excited energies (even with the aid of the variational 
method), we need to use finer temporal lattice spacings. This allows us to use as many time slices 
of information as possible to reconstruct more accurately the signal due to a particular state, since 
the signal in the Euclidean space declines exponentially as Et. In principle, we could use very fine 
lattices while keeping the volume big enough to avoid the "squeezing" systematic effect; however, 
the computational cost scales significantly (power of 5-6) with the lattice discretization. We adopt 
the approach of anisotropic clover lattices to keep the spatial lattice spacing coarse, avoiding finite- 
volume systematic error, and to make the temporal lattice spacing fine enough to extract towers of 
nucleon excited states. However, tuning the 0(a) -improved parameters for fermion actions in the 
dynamical gauge generation correctly is a much more difficult task for anisotropic lattices. We use 
Symanzik-improved gauge action and clover fermion action with 3-dimensional stout-link smeared 
gauge fields; the gauge ensemble is generated using the (R)HMC algorithm. Our spatial lattice 
spacing is a s = 0.1227(8) fm (determined using mjj), and the renormalized anisotropy t, R = a s /a t 
is 3.5. Ref. ^ shows a detailed study of the dynamical anisotropic lattice parameter settings, and 
Ref. [Q] reports basic lattice properties along with the ground-state hadron spectrum. 

Since then, we have also moved on to generate m n = 230 MeV 24 3 and 32 3 lattices. A stream 
of m n = 180 MeV ensembles is slowly progressing (for both 24 3 and 32 3 volumes). In Ref. [Q], we 
demonstrated and compared various strange-quark mass-setting approaches. We tuned the param- 
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Figure 2: (Left panel) The location of the dynamical ensembles used in this work in the sq-Iq. plane. 
The leftmost circle (black) indicates the physical point {l^ ys , s a hys }, while the red, green and blue points 
are the lattices with a t m s = —0.0540,-0.0618,-0.0743 respectively in Ref. [H], The horizontal dashed 
(pink) line indicates constant sq at the physical point, and the diagonal line indicates three-flavor degenerate 
theories. The purple points are from recent measurements of the 230 Me V ensemble, showing that it does not 
deviate away from the Nf — 3 value. Mass-ratio chiral extrapolations as functions of the h and sz for octets 
(center panel) and and sq for decuplets (right panel). The lines indicate the "projected" leading chiral 
extrapolation fits in h and Ss (Iq and sq) while keeping the other one fixed. The black (circular) point is the 
extrapolated point at physical h and s? (la and sq). The masses for N, £ and A are 0.962(29), 1.203(11), 
1.122(16) GeV, respectively, and the masses for A, E* and S* are 1.275(38), 1.426(25) and 1.580(13) GeV, 
respectively. 



eter sq = 4*7 as close as possible to the corresponding experimental value at the SU(3)f- 
symmetric point (in our case, it is the 875-MeV-pion ensemble). When we reduce the sea-quark 
mass, the sq parameter remains roughly constant as low as the 383 MeV pion ensemble. (Fur- 
thermore, the sq parameter is more sensitive to the strange-quark mass than other methods, such 
as /-parameter. Since it is a ratio, there is no need to worry about estimating the shift in the lat- 
tice spacing when decreasing the sea-quark mass or generating more statistics.) We measure the 
same quantity on the 230-MeV 24 3 ensemble and sq does not deviate from the expected value. 
The sq parameter is a stable and useful observable for setting the strange-quark mass; we have not 
observed any notable deviation for a wide range of sea-quark masses. 

In Ref. [||], we reported a naive extrapolation in mass ratios (using vyiq or m; as the reference 
mass) to obtain the meson and baryon masses. Note that at each sea-strange mass, we only use the 
larger-volume set for the lightest pion mass to avoid large finite-volume effects. The advantage of 
using the mass ratios instead of the masses is that we gain smaller statistical (and possibly system- 
atic) errors than the mass itself due to cancellation of systematics and removal of the ambiguity 
associated with setting the lattice spacing. Using the same data, we update the mass-ratio chiral 
extrapolations by modifying the next-to-leading-order heavy-baryon chiral perturbation theory. We 
find the finite-volume corrections are negligible according to ChPT estimates and the extrapolated 
baryon masses agree with experimental values significantly better than with a linear extrapolation. 
Preliminary results are shown in Fig. 

3. Methodology: Cubic-Group Operators and Distillation 

The Hadron Spectrum Collaboration (HSC) has been investigating interpolating operators pro- 
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jected into irreducible representations (irreps) of the cubic group [|7|] in order to better calculate 
two-point correlators for nucleon spectroscopy. Baryon correlation functions were evaluated using 
the displaced-quark operators described in Refs. [[7|] and employed in spectrum studies in Refs. [j|]. 
In the cubic group Oh, for baryons, there are four two-dimensional irreps G\ g ,G\ u ,G2 g , G% u and 
two four-dimensional irreps H g and H u . (The subscripts "g" and "u" indicate positive and nega- 
tive parity, respectively.) Each lattice irrep contains parts of many continuum states. The Gi irrep 

I "7 Q 1 1 ^ S *7 Q 

contains J = ~ , £ , | , -5-, . . . states, the H irrep contains J = i, |, £ , |, . . . states, and the G2 irrep 
contains J = |, 2, ■ ■ ■ states. The continuum-limit spins / of lattice states must be deduced by 
examining degeneracy patterns among the different Oh irreps. Using these operators, we construct 
an r x r correlator matrix and extract individual excited-state energies by fitting with single- and 
multiple-exponential functions. 

However, the calculations using these cubic-group operators require multiple orientations in 
order to maximally overlap with a wide range of quantum numbers, which is quite expensive. 
Furthermore, as we go to lighter and lighter pion masses, there will be more decay modes open, 
even for the lowest energy at a specific quantum number. We need to extend the matrix to include 
the multiple-particle operators (so that we can further understand the nature of the "state" in our 
calculation) and "disconnected" operators. Further, we need to achieve better precision for each 
state to distinguish among them. 

A new way to calculate timeslice-to-all propagators, "distillation", has been proposed in Ref. [^]. 
The method is useful for creating complex operators, such as those used in the variational method, 
allows the operators to be decided after performing the Dirac inversions and reduces the amount 
of time needed for contractions. Distillation uses color-eigenvector sources to improve on noisy 
estimators, giving better coverage of relevant degrees of freedom. Increasing the number of "hits" 
improves statistics faster than 1/yW. The method can combine with stochastic methods, which 
might be desirable if the number of sources needed to cover the volume becomes too large. 

The distillation operator on time-slice t can be written as 

U(t)=V{t)V\t) -» LVO = £ vP(t)4 k) Ht), (3-D 

k=\ 

where the V(t) is a matrix containing the first through k th eigenvectors of the lattice spatial Lapla- 
cian. The baryon operators involve displacements (£^,) as well as coefficients (S ai0 , 2a3 ) in spin 
space: 

where the color indices of the quark fields acted upon by the displacement operators are contracted 
with the antisymmetric tensor, and sum over spin indices. Then one can construct the two-point 
correlator; for example, in the case of proton, 

4 2) [z d , x u , z u ] (t>, t) = <&('■■;■■*) (t')rP (t', t)d ^(t>Mu*\t>,t)^>(t) 

-^V)^^ (3.3) 

where the "baryon elemental" 

= ^ (^iv (,) ) a (®2V (i) ) b S aiaia3 (3.4) 
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can be used for all flavors of baryon and quark masses with the same displacements on the same 
ensemble, and the "perambulator" 



[t',t)=vHt')M a l p (t',t)V(t) (3.5) 



can be reused for different baryon (and meson) operators after a single inversion of the M matrix. 
There is a large factor of computational power saved by "factorizing" correlators in terms of el- 
ementals and perambulators, and we can use the same elementals and perambulators to contract 
various different correlators. 

The distillation method is demonstrated on four nucleon G\ g operators, two of which are local 
to a single site and the remaining two having a singly displaced quark field. The corresponding 
4x4 matrix of correlators was computed on 316 configurations of the 16 3 x 128 (« 380 MeV) 
lattice ensemble. Modeling the correlator noise-to-signal ratio with a + as a function of the 



number of distillation eigenvectors included (N in Eq. |3.1| ) gives a best-fit exponent p ~ 1 . 1 (2); that 
is, increasing the number of vectors decreases the noise considerably faster than simple statistical 
scaling. The variational method [ |l0| ] is used to extract the masses of the lowest two states in the G\ g 
spectrum. The extracted mass dependence on the number of eigenvectors is shown in the left panel 
of Fig. ||. We find consistent masses among the larger numbers of eigenvectors, with an increase in 
statistical precision as the number of vectors is increased. 

A procedure called "pruning" is used to reduce the number of operators included. A practical 
procedure is to calculate the correlators with the same source and sink operator (i.e. the diagonal 
elements of the full correlator matrix) and sort them according to their "individuality". One way to 
systematically prune is to take the matrix of inner products of the effective masses of each correlator 
across all time slices and sort them from there. The middle panel of Fig. || shows a subset of various 
Gig displacement-operator correlators sorted by the values of their inner products with respect to 
one another. Any overlap larger than 70% is marked by yellow, while the remaining values are 
depicted as magenta to blue colors. (G\ g has the largest overlap amongst its operators of all the 
baryon irreps.) Notice that the inner-product matrix has a clear block structure when organized this 
way. Each yellow block consists of a set of operators that yield nearly identical effective masses 
as functions of time. We exclude many operators by such a selection process and pick a total of 24 
operators among different quark orientations for each irrep. Then we refine to a smaller sub-matrix 
of 8 x 8 by filtering the matrix by condition numbers; we perform a variational-method analysis on 
the resulting 8x8 matrix and retain the lowest 4 eigenstates. 

Fig. |3| shows preliminary results for nucleon spectroscopy with pion mass around 380 MeV 
on a 16 3 x 128 volume, using "distillation" []|] with N = 32 eigenvectors. We observe a similar 
distribution of states as the previous study. Further measurements on the larger volume and inves- 
tigation of decay thresholds and potential two-particle states are underway. (See the proceeding by 
K. Juge, reporting two-particle results using distillation.) 



4. Conclusion and Outlook 

We report an ongoing effort to solve the mysteries of baryon resonances. The ground-state 
baryon masses are in reasonable agreement with experiment and consistent among different groups 
with different actions if systematics are taken into consideration. Our new technique, "distillation", 
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Figure 3: (Left) The fitted energies for the first-excited (top) and ground state (bottom) as functions of the 
number of distillation eigenvectors N. (Middle) A subset of G\ g irrep operator correlators, grouped by their 
inner product. The yellow blocks indicate overlap > 70%. G\ g is the worst case among all the irreps. (Right) 
Nucleon excited spectrum sorted according to cubic-group irrep. 



will greatly improve precision in our future calculations for extracting excited-state masses using 
cubic group-irrep operators, which provide powerful probes to extract highly excited resonances. 
A preliminary result for nucleons on Nf = 2 + 1 m % = 380 MeV is shown in this proceeding. Work 
on larger volumes (with a modified stochastic distillation) are under development. Meanwhile, 
parallel work from the HSC for meson spectroscopy with exotic quantum numbers and baryons 
using derivative operators are also in progress. Multi-particle operators are under investigation 
to distinguish these from resonances. We are also investigating the application of the distillation 
method to form factors to help us understand the nature of specific states. 
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